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Abstract 

Using a continuum bead-spring Monte Carlo model, we study the anomalous diffusion dynamics 
of a self-avoiding tethered membrane by means of extensive computer simulations. We focus on the 
subdiffusive stochastic motion of the membrane's central node in the regime of flat membranes at 
temperatures above the membrane folding transition. While at times, larger than the characteristic 
membrane relaxation time rpt, the mean-square displacement of the center of mass of the sheet, 
(Re), as well as that of its central node, (i?^), show the normal Rouse diffusive behavior with a 
diffusion coefficient scaling as Djy (x N^^ with respect to the number of segments N in the 
membrane, for short times i < we observe a multiscale dynamics of the central node, {Rn) oc t", 
where the anomalous diffusion exponent a changes from a ~ 0.86 to a ~ 0.27, and then to a ~ 0.5, 
before diffusion turns eventually to normal. By means of simple scaling arguments we show that 
our main result, a ~ 0.27, can be related to particular mechanisms of membrane dynamics which 
involve different groups of segments in the membrane sheet. A comparative study involving also 
linear polymers demonstrates that the diffusion coefficient of self-avoiding tethered membranes, 
containing segments, is three times smaller than that of linear polymer chains with the same 
number of segments. 

PACS numbers: 87.16.D-, 02.50.Ey, 87.15.A-, 87.15.Vv 
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I. INTRODUCTION 



There has been considerable interest recently in understanding the statistical properties of 
polymerized (or tethered) membranes [l| . This interest in large part is due to the membrane 
behavior which is much richer than that of polymers, their one- dimensional analog. In addi- 
tion, this interest is justified by a variety of real systems like red-blood-cell cytoskeletons Q], 
graphite oxide sheets js, 4] or dispersed silicate (clay) platelets [s], 6| which can be modeled by 
networks of fixed connectivity, generally referred to as polymerized membranes. Along with 
the experimental studies, self-avoiding polymerized membranes have also attracted remark- 
able interest from the point of view of basic research in recent years. Their static properties 



10 



11 



12 



3,Q, 



15 



16 



17, a 



have been studied analytically and numerically j?, s], 0, 
Much of these studies have been spent in the pursuit of the so called "crumpling transition" 
between a low-temperature fiat phase and a high-temperature crumpled phase until it was 



realized 



2l| that self-avoiding membranes are always flat (with an infinite persistence 



length), i.e., their radius of gyration Rg scales with linear size L a.s Rg cc L'^ where the Flory 
exponent u ^ 1. The flat phase arises even without explicit bending rigidity because the 
resistance to in-plane shear deformations leads to anomalous stiffening of the surface in the 
presence of thermal fluctuations. 

In contrast to static properties, the membrane dynamics is less well-understood. Earlier 
analytical and numeric studies 0, Q, 



22 



23[ | have revealed that the self-avoiding restric- 



tions considerably modify the relaxation times of the tethered surface. Thus the typical 
relaxation time r^j of a tethered membrane in the case of Rouse dynamics when hydrody- 
namic interactions are neglected has been predicted by simple scaling arguments 7] to vary 
as Tji oc OC i?f If, as in polymer physics, one introduces a dynamic exponent z, 

describing the relaxation process as r^j oc i?^, then one gets z = 2 + 2/v (for linear polymers 
one has z = 2 + 1/z/) . Usually, is considered to be the time needed for the membrane 
to diffuse its radius of gyration. For tethered membranes, highly permeable to the solvent 
as in isolated spectrin networks, one is in the Rouse regime j22| and the diffusion coefficient 
Z^AT scales with membrane size L as oc L^^. Thus the time it takes for such a flat mem- 
brane to move a distance Rg is proportional to L^. In contrast, for impermeable membranes 
(like, e.g., erythrocytes) where solvent backflow (i.e., a long-ranged hydrodynamic interac- 



tion) is important, one has in d-dimensions in the case of Zimm dynamics th oc R'^ oc L 



dv 
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I.e. z 



d) and oc L ^. Thus one may view permeability as constituting two differ- 



ent dynamic universality classes of tethered membranes [22] whereby these classes (Rouse 
dynamics - highly permeable membranes, or Zimm dynamics - impermeable membranes) 
are observed for a wave vector independent (or, dependent) friction coefficient. Recently, a 



series of simulation studies by Pandey et al. [2j, |25|] has revealed a multiscale stochastic 
dynamics of tethered membranes at times before normal diffusive behavior is reached. The 
displacement motion of the central node, Rn, of a four-coordinated coarse-grained model 
membrane has been observed to undergo a subdiffusive mean-square displacement (MSQD) 
(Rn) oc with the exponent a attaining different values in the short and intermediate time 
regimes before turning to normal diffusion for t > tji with a = 1. 

In the present work we employ an efficient off-lattice Monte Carlo algorithm, focusing on 
the subdiffusive dynamics of self-avoiding tethered membranes and comparing some of the 
salient dynamic features to those of linear polymers. Our observations, based on extensive 
computer simulations, largely confirm those of previous investigators 2J, |25l]. As a step 
forward, however, we suggest a scaling theory which explains our findings for the anomalous 
membrane dynamics, relating the observed values of a to the specific stochastic motion of 
particular groups of sheet segments. 

After briefly introducing our model in Section II, this is considered in Section III where 
we focus on the main results of our investigation and their interpretation. We close this 
paper with a brief summary and discussion in Section IV. 



II. MODEL SYSTEM AND SIMULATION PROCEDURE 

We study a coarse-grained model of self- avoiding tethered membranes, embedded in three- 
dimensional space. The membranes have a hexagonal lattice structure where each monomer 
interacts with six nearest-neighbors - Fig. [H There are altogether = (3L^ — 3L + 1) 
monomers in such a membrane where by L we denote the number of monomers on the edge 
of the network (i.e. L is the linear size of the membrane). In this model, spherical particles 
of diameter a are connected in a fixed geometry by flexible strings of length /. To prevent 
self-intersection of the membrane, the maximum length of the strings between the centers of 
the spheres must be / < VSa, then the membrane is self-avoiding in that it cannot intersect 
itself. 
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FIG. 1: (Color online) A snapshot of a tethered membrane at T = 1.0 with linear size (the edge 
length of a regular hexagonal sheet) L = 50 which contains N = 7351 monomers. Periphery 
segments at the rim of the membrane are shaded grey (green) while the six monomers at the 
vertices (corners) of the sheet are dark grey (red). 

The bonded nearest neighbor monomers on the membrane interact with each other 



through the finitely extensible nonlinear elastic (FENE) potential 26| where a bond / has a 
maximum length /max and a minimum length /min, 



UfENe{1) = --^(^max - ^o) 



1 



l-lo 

^max ^0 



(1) 



The minimum of this potential occurs for I = Iq , ?7fene(/o) = 0, near /q it is harmonic, 
with K being a spring constant, and the potential diverges to infinity both when / /max 
and when / /min- Choosing our length unit /^ax = 1-0, we choose the other parameters 
as /min = 0.2, /q = (/min + ^max)/2 = 0.6, K/ksT = 5, where T denotes the absolute 
temperature, and ks is the Boltzmann constant. 

Self-avoidance is observed by the interaction between particles which are not nearest 
neighbors on the network. The nonbonded interaction between monomers is described by a 
Morse potential where r is the distance between the monomers, 

^^Morse(?^)/eM = exp[-2a(r - rmin)] - 2exp[-a(r - rmin)] (2) 

with parameters eM/ksT = 1, a = 24. The minimum of this potential occurs for r = 
'^min, ^Morse('^min)/eM = "l- For a = 24, UMorscir) essentially is zero for r > 1.25 rmm- 
Choosing then units of length such that rmin = 0.8, we hence can take t/Morscl'^ > 1) = 0. 
The repulsive part of this potential guarantees self-avoidance of the membrane. 

We have used the standard Monte Carlo procedure to investigate the thermodynamic 
properties of self-avoiding tethered membranes. The total energy (Hamiltonian) is the 
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sum of Eqs. ([T]) and ([2]). In each Monte Carlo update, a monomer is chosen at random 
and one attempts to displace it randomly by displacements chosen uniformly from the in- 
tervals —0.25 < Ax, Ay, Az < +0.25. The attempted move is accepted or rejected ac- 
cording to the conventional Metropolis criterion by comparing the transition probability 
W = exp{—AU/kBT) (where AU is the energy difference between the configurations after 
and before the trial move) with a random number uniformly distributed between zero and 
unity. If W exceeds this random number, the attempted move is accepted, otherwise it is 
rejected. Time is measured in Monte Carlo steps (MCS) per monomer whereby a single MCS 
is elapsed after monomers are picked at random and given the chance to perform a trial 
move. Since our potentials are constructed such that the membrane cannot intersect itself 
in the course of random displacements of monomers, one does not need to check separately 
for entanglement restrictions during the simulation. Thus the algorithm is reasonably fast. 
Nevertheless, the simulation takes quite a long time for large self-avoiding membranes to 
equilibrate and then move a substantial distance in space. This and the necessity to attain 
very good statistical accuracy have limited our investigations to sizes L < 50. Eventually, 
we would like to note that the interactions used in the present off-lattice model, albeit some- 
what more refined and complicated than the simple potential used in earlier simulations on 
a cubic lattice 



25l |. do not change the physics of the problem and lead qualitatively to 



the same results. 



III. RESULTS 

Before we focus on the subdiffusive dynamics of our membranes, we show in Fig. [2] the 
scaling behavior of the gyration radius, 

1 ^ 

i=l 

N 

where is the position of the ith monomer of the membrane while rem = Ji ^^^i its 

i=l 

center of mass location. We also sample the eigenvalues A^^^, A^g^j, X^^^ of the inertial 
tensor, 

1 ^ 

Ic^P = j;^T.(^^-^cJ{rf-rL), (4) 
1=1 
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FIG. 2: (Color online) Log-log plot of the mean-squared radius of gyration, i?^, and the three 
eigenvalues, A^a,x' -^med' ^min against the linear size L of flat tethered membranes (5 < L < 50) 
at r = 1.0. The respective exponents are given in the legend. In the inset we show the variation 
of membrane asphericity, A = A^j^/A^a^^ which vanishes with L as A (x L^^-^^. 
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where G {x,y,z}, the sum is taken over all particles of a given configuration, and r[ 
is the a component of the center of mass radius vector for a given configuration. The three 
eigenvalues are ordered according to magnitude A^^^^ > A^^^^j > A^j^. The directions of the 
principal axes are given by the three eigenvectors corresponding to the three eigenvalues. 
For a planar membrane, the eigenvector associated with A^j^ is perpendicular to the plane 
of the membrane while the eigenvectors associated with A^^^^^ and A^^^ lie in the plane of 
the membrane. It is evident from Fig. [2] that our membranes are indeed fiat with scaling 
exponents u = L02±0.01 for Rg, v = L02±0.01 for A^^^, and u = L06±0.01 for A^.^- The 
asphericity ratio A = (A^jn/A^j^^) oc L"'^'^^ tends to zero with an exponent ua = 0.32 ±0.01, 
indicating that these membranes are indeed asymptotically fiat. Note that these data have 
been obtained at T = 1.0 well above the temperature of the first folding transition j2^, 
= 0.89 ± 0.01. In Fig. [2] and in the following figures the error bars do not exceed the size 
of the symbols. 

Turning now to membrane dynamics in the Rouse regime, one may assume that each 
segment of the membrane moves under the influence of surface forces (surface stretching due 
to near-neighbors and excluded volume forces due to distant neighbors), and a random force 
representing thermal noise. As far as the contribution of inertial terms to membrane motion 
can be neglected for sufficiently long times, one may assume that the relevant djTiamics is 
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purely diffusive. With p being the rate of position changes of monomers per unit time and 
z = 2/z/ + 2, the dynamic exponent, one may write the relaxation time of the membrane tr 
as 

TR = p-'Rl = p-'N'^/' (5) 

If monomeric orientations add up randomly and one neglects correlations, the MSQD of the 
membrane center of mass is 

93{t) = {[fcUt) - rUO)?) = P{ (-(^) )Nt = p^t, (6) 

because each monomeric motion moves the center of mass by a random displacement of 
the order l/N, I being the bond length. There are pN such random motions per unit time. 
Invoking the Einstein relation gsit) = 2dD]^t (where d is the spatial dimensionality), one 
thus concludes 

DN(xp{f)/N. (7) 
From Fig. [3] it is evident that this prediction, Eq. ([7]), is indeed nicely confirmed by the 
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FIG. 3: (Color online) Variation of the diffusion coefficient Dj\f of a membrane with the number 
of monomers for membrane sizes 5 < L < 50 at T = 1.0. The measured slope is —1.01 it 0.02. 



simulation. The relaxation time tr and the scaling law 2;z/ = 2z/ + 2 is then understood by 
the condition that the membrane is relaxed when its center of mass has diffused over its own 
size Rg, i.e.. 



gsiTR) oc p^-^TR oc (f )iV-/2-i oc i?J oc {f)N'' 



(8) 



whence zu = 2 + 2^ follows. 

These arguments can be carried over for the MSQD of the membrane central node as 
well. We define 



g^{t) = {[f^{t)-r^{0)r)cx{P){ptr (9) 

at time t < tr anticipating that the central node exhibits anomalous diffusion with an 
exponent a < 1. For short times (pt) < 1, of course, a nearly free diffusion of the central 
node takes place, and thus gi{t) for pt ^ 1 should be of the order of oc (P). Requiring now 
that gi{Tji) ^ _R^, one gets a scaling relation for a, 

giim) oc {P){pTRr oc (f)iV"-/2 oc {f)Nr (10) 

Thus for flat membranes with z/ = 1, one has = 1 + z/~^ = 2, and one would then expect 
to see a time interval t < where gi{t) oc t^/^. Additional information for the subdifFusive 
dynamics of the membrane may be obtained if one defines in analogy with the case of linear 



polymers 



28l | the MSQD of a central node measured in the center of mass coordinate system 



of the membrane, 

g2{t) = {[rn{t) - rcUt) - UO) + rcmm^), (11) 

and also for the averaged MSQD of the six monomers at the vertices of the hexagonal sheet 
in the laboratory system of coordinates, 

9^it) = {li2l'-'^it)-m?)■ (12) 

i=l 



In the center of mass coordinate system of the membrane, 

6 

'6 



1 ^ 

g5{t) = (- J2m) - r,m{t) - f,(0) + rcmiO)?)- (13) 



Evidently, for t < tr one should observe g2{t) ~ giif) and g^it) ^ giit) whereas g2{t) oc -Rg 
for t ^ Tji since the central monomer cannot travel farther from the center of mass than 
the membrane size, of course. 

A general impression about the time variation of the various MSQD gi{t) and the simi- 
larity in the stochastic dynamics of linear polymers and tethered membranes may be gained 
from Fig. |H It is seen that the course of the functions gi with time is qualitatively very 
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FIG. 4: (Color online) Log-log plot of mean square displacements gi{t), g2{t), g3{t), Qiii), g^{t) at 
T = 1.0 plotted vs time t (measured in Monte Carlo steps) for a linear polymer of length N = 256 
(a), and for a membrane of size L = 10 which contains N = 271 monomers (b). Dashed lines 
indicate the scaling behavior of the central segment MSQD, gi{t) oc t", and of the center of mass, 
33 (t) oc t, with elapsed time while horizontal lines denote the time averages of the radius of gyration 
Bj^ and of the "end-to-end" distance of the chain, (a), and of the membrane, i?gg (b), the latter 
being measured as the distance between the opposite vertices of the hexagonal sheet. Evidently, 
cf. (b), by defining r/j from g^ijii) = as the mean relaxation time of the membrane one can 
verify that gi{TR) ^ gsirn). 

similar for both linear polymers and tethered membranes. In both cases, cf. Figs. 111(a) 
andlD^b), one finds that the center of mass performs normal diffusion with gsit) = QD^t. 
One can, therefore, compare the relative diffusivity of polymers and membranes, contain- 
ing the same number of monomers A^, say a chain with N = 256 and a membrane with 
L = 10, i.e., = 271, with identical forces acting between the repeating units. Our anal- 
ysis shows that in a good solvent, T = 1.0, one obtains GD^v (polymer) = 6.6 x 10~^ and 
6/^ AT (membrane) = 2.2 x 10^^, i.e., a linear self-avoiding chain of segments moves three 
times faster than a self-avoiding flat sheet in the case of Rouse dynamics. We find this re- 
sult rather remarkable since this decrease in mobility is solely and entirely due to the higher 
topological dimensionality of the membrane. 

A marked difference between chains and membranes, however, is revealed if one looks at 
the subdiffusive behavior of the central monomer in both cases. For times shorter than the 
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typical relaxation time, t < tr, the central node of the polymer chain is observed to diffuse 
like gi{t) oc t^'^^ (i.e. very close to the expected t^'^^ power law) while for the membrane 
one finds a much smaller power gi{t) oc t^'^, seen also by Pandey et al. 2^- It might be 
argued that this small exponent a ^ 0.3, describing the subdiffusive behavior of a tethered 
membrane, reflects a membrane-specific dynamic mechanism which shows up at t < tr. In 
the following we suggest a possible interpretation and a simple scaling derivation for the 
observed value of this novel exponent a. 



g,(t) 
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FIG. 5: (Color online) Log-log plot of the central node MSQD 171 vs time t at T = 1.0 for a 
membranes with linear size 5 < L < 50. All data are averaged over 100 simulation runs. Dashed 
lines denote power law variation with different exponents (see legend) corresponding to the various 
subdiffusive regimes. 

We first look more closely at the main data of our study showing the MSQD of the central 
monomer of a tethered membrane in Fig. O The different regimes of subdiffusive motion 
of the central node are indicated by power laws with exponents, specified in the legend 
of Fig. [5l It is seen that for very short time, < t < 1, each segment indeed performs 
displacements which are not constrained by the topological connectivity of the network so 
with gi{t) oc t^'^^ one observes an extremely shortlived nearly normal diffusion. At late 
times, t > Tji, the normal diffusive motion sets on eventually, and gi(t) ^ g^it) oc GD^t. We 
should like to point out that at late times the averaging of the correlation functions gi for 
the two largest system sizes, L = 30 ^ 50, is not perfect due to a progressively deteriorating 
statistics, however, it is beyond doubt that their ultimate slope corresponding to normal 
diffusion should be unity. In the intermediate time interval our data yields a subdiffusive 
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motion of the central monomer with gi{t) oc where a = 0.27 ± 0.01. Due to strong 
finite-size effects this value of the a can be unambiguously established for sufficiently large, 
20 < L < 50 membranes only. We note that very close values for a ~ 0.25 -r- 0.32 have been 
observed recently in the computer experiments of Pandey et al. {25 1 for the case of tethered 
membranes in the good solvent regime at temperatures 2.0 < T < 10. Such behavior cannot 
be explained by means of the exponent a = 0.5 which follows from the estimate, Eq. (fTOl) . 

We believe that a possible explanation of this sluggishness of flat membranes at early times 
1 < t < Ti < T/j may be found if one assumes that in this interval only the most loosely 
bound monomers (those at the membrane periphery, or rather, those at the six vertices of 
the hexagonal sheet) actually contribute to a displacement of the membrane center of mass 
while all monomers with six-fold coordination in the bulk of the membrane are virtually 
blocked by their neighbors and for this short time hardly move. As far as the membrane 
retains its flat shape and does not fold, the maximal displacement of these loosely bound 
monomers cannot exceed the effective thickness Amin- Indeed, a comparison of Fig. [2] and 
Fig. [5] shows that the MSQD, performed by a membrane of linear size L during the time ti 
(the latter is given by the intersection point of the tangent y oc to gi, and gs oc t) amounts 
to (73 (ti) ~ A^jjj. One can, therefore, estimate the characteristic time ti if, in analogy to 
Eq. ([6]), one considers 

N 



9-^(^) = Pi{i] (14) 



in case that only the membrane periphery of length oc y N contributes to the center of mass 
displacement. Thus during 1 < t < ti the diffusion coefficient of the membrane becomes 
D oc N"^/"^ oc reflecting the slow displacement of the center of mass. Eq. (HM therefore 
suggests Ti oc L^^^'^™'" ^ L^-^. With 5'i(ti) oc L"^^^^'^""") ^ A^j^ one obtains then the broken 
exponent for a periphery- driven membrane ap = (2z/jnin)/(3 + 2z/jnin) = 0.32. 

In contrast, if only a flnite number of loosely bound monomers at the vertices effect the 
net displacement of the center of mass, one obtains then 

and therefore temporarily D oc A^~^ oc so that ri oc /,4+2iyniin ~ L^'^. In this case one 
gets the exponent for anomalous diffusion of a vertices- driven membrane ay = (2z/min) /(4 + 
2u^^) = 0.26. 
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Thus we obtain two estimates which may be considered as the lower and upper bounds 
of the anomalous diffusion exponent, < a < ap, depending on the particular mechanism 
involved in the diffusive motion. The measured value of a ~ 0.27 lies indeed within these 
limits. Of course, one should bear in mind that most probably neither mechanism of diffu- 
sion (vertices-driven, or periphery-driven) takes place alone and the real process involves a 
mixture of both. Moreover, at this point we cannot rule out the possibility that at times 
T~i ^ t < tr all membrane segments eventually get the chance to perform an elementary 
move and thus contribute to the center of mass motion. Such a possibility would imply that 
during this time interval of subdiffusive motion one observes a MSQD gi{t) oc with an 
exponent a = 0.5, cf. Eq. fllOp . We have indicated such a behavior in Fig. Elby a dashed 
line with slope 0.5 and it appears compatible with the course of gi{t) in between ti and 
for our largest membranes L = 30 ^ 50. If such a diffusive regime really exists, it would 
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25| . It is clear, however. 



underline the multiscaling character of tethered membranes 
that larger micelles need to be simulated with satisfactory statistics before an unambiguous 
conclusion in this respect can be drawn. 



IV. SUMMARY AND CONCLUSIONS 



In the present work we have studied the stochastic dynamics of flat self-avoiding tethered 
membranes which are assumed to be completely permeable to the surrounding good solvent 
and are thus expected to display typical Rouse behavior. By means of extensive Monte Carlo 
simulations we find that the static properties of our tethered membranes are described by 
scaling exponents which agree very well with the appropriate theoretically predicted values. 
Thus the radius of gyration scales with membrane linear size L as i?^ oc L"^" with u = 1.02 ± 
0.01, and the membrane thickness, A^j^ oc L^'^™'" with roughness exponent z^min = 0.70 ±0.01 
while the membrane asphericity vanishes asymptotically as A = A^j^/A^j^x ^ L^'^^-'^ with 
ua = 0.32 ±0.01. 

In the regime of Rouse diffusion we find with good accuracy that the diffusion coefficient 
(X N~'^, as predicted, whereas the typical relaxation time of such polymerized membranes 
grows as r/j oc with the linear dimension L. A comparative study, involving linear 
polymers too, reveals also that the diffusion coefficient of permeable self-avoiding tethered 
membranes, containing segments, is three times smaller than that of linear polymer chains 



12 



with the same number of segments. 

Our main concern in this study, however, is with the subdiffusive motion of the membranes 
central segment at times t < t/j. Our numeric studies reveal several regimes of anomalous 
diffusion whereby the central node MSQD grows as gi oc t^'^^ for t < 1, then gi oc t'^'^^ 
for 1 < t < Ti, further, with gi oc t^'^ at time Ti < t < r^, before turning eventually to 
normal diffusion with gi oc t for t > Tn. We use simple scaling arguments to interpret our 
observation and suggest that the anomalous diffusion exponent a ~ 0.27 which we find in 



agreement with recent studies [2J, |25(] most probably reflects several particular mechanisms 
of membrane motion. These mechanisms involve different groups of loosely bound membrane 
monomers whose random hops predominantly contribute to the center of mass motion of 
the whole membrane at times when most of the inner monomers are mutually blocked by 
their nearest neighbors and, therefore, remain rather immobile. The particular geometry of 
the membrane sheet (e.g., square, hexagonal, or rhombic) is expected to enhance the role 
of either periphery, or vertex monomers, and therefore slightly modify the observed value of 
the anomalous exponent a according to Eqs. iHM and ffTSl). T his would explain some small 
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25|. 



deviations of our data from that of earlier measurements 

We believe that our results shed some light and provide insight into the complex dynamics 
of polymerized membranes. It is, however, clear that further work is needed before the nature 
of the membrane stochastic dynamics is definitely established and understood. 
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